Two-point correlation properties of stochastic "cloud processes" 
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We study how the two-point density correlation properties of a point particle distribution are 
modified when each particle is divided, by a stochastic process, into an equal number of identical 
"daughter" particles. We consider generically that there may be non-trivial correlations in the 
displacement fields describing the positions of the different daughters of the same "mother" particle, 
and then treat separately the cases in which there are, or are not, correlations also between the 
displacements of daughters belonging to different mothers. For both cases exact formulae are derived 
relating the structure factor (power spectrum) of the daughter distribution to that of the mother. 
These results can be considered as a generalization of the analogous equations obtained in .1\ for 
the case of stochastic displacement fields applied to particle distributions. An application of the 
present results is that they give explicit algorithms for generating, starting from regular lattice 
arrays, stochastic particle distributions with an arbitrarily high degree of large-scale uniformity. 
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I. INTRODUCTION 

Point processes, i.e., stochastic spatial distributions of 
identical point particles, provide a very useful mathe- 
matical scheme for many different A'^— body and granu- 
lar physical systems such as crystals (both regular, and 
perturbed) @, SSl, quasi-crystals 0, structural glasses, 
fluids [y| , self-gravitatingsystems in astrophysics and cos- 
mology (see, e.g., 0, SQ)- They find also many appli- 
cations in a wide range of scientific fields: computer im- 
age problems [13], and bio-metrical studies ill,] are only 
some examples of systems which are usually represented 
as specific point processes with appropriate spatial cor- 
relation properties. The extension of knowledge about 
this class of stochastic processes can therefore be of fun- 
damental importance for the description and analysis of 
many scientific topics. Indeed this is the reason why con- 
siderable mathematical effort has already been invested 
in the study of this class of systems, and many useful 
results have been derived (e.g., see [l3. [l3. [l4| ) . 

In this paper we present the equations for the structure 
factor (or power spectrum) of a point process obtained 
as follows. We start with an arbitrary uniform spatial 
point particle distribution with a known structure factor 
(SF). We now suppose that each of these "mother" parti- 
cle splits into a cloud of m identical "daughter" particles, 
where to is a cloud-independent constant. Each daughter 
particle is then assumed to be displaced from its mother 
position by a stochastic displacement which may, or may 
not, be correlated with the displacements of other parti- 
cles. In other words each set of m particles initially lying 
at the same spatial point "explodes" forming a "cloud" 
of particles around it. For this reason we call the point 
process so generated a cloud process. 

We suppose that the displacements applied to differ- 
ent particles belonging to the same mother are symmet- 



rically distributed with arbitrary pair correlations. One 
can choose, for instance, these correlations in order to fix 
certain moments of the mass dispersion of each cloud. We 
will distinguish in the following between the two cases in 
which the displacements applied to particles originating 
from different mothers are, or are not, correlated. We 
note that the results obtained for these cases are gener- 
alizations of those obtained in [11], where the case of a 
single daughter for each mother was analyzed and solved 
exactly. 

The main immediate application we have in mind of 
the present study, and the one we discuss at some length, 
is the determination of the constraints on the stochastic 
displacement field and the number of daughters m re- 
quired to produce a target behavior of the SF at large 
scale (i.e. small wave number k) in the daughter particle 
distribution. More specifically we are interested in the 
case of a large scale SF inherent to superhomogeneous 
(or hyper-uniform) stochastic point particle distributions 
(iSl [la . UJ^ . This class of distribution is characterized 
by the convergence of the SF to zero as fc ^ 0. We will 
show explicitly that, for the case that the mother distri- 
bution is a regular lattice array, the associated (positive) 
exponent in the fc ^- limit of the SF is related to the 
conservation of the local mass moments in the passage 
from the single point particle to the cloud of daughter 
particles^. 

The paper is structured as follows. In the next section 



^ We note that approximate heuristic derivations of some of these 
results can be found in the cosmological literature (see, e.g., 
[Si, lg| ) in which constructions of this kind are considered in dis- 
cussions of "causality bounds", i.e., bounds imposed on the large 
scale behavior of the SF by the existence of a causal horizon in 
cosmological models. 



we introduce our notation and the essential definitions, 
and then derive a general expression for the SF of the 
daughter distribution, before averaging over the realiza- 
tions of the cloud process. In Sect. IIIII we derive our 
general result for the SF for the case that there are only 
correlations between displacements of particles derived 
from the same mother. In the following section we apply 
this result to the specific case that the mother distribu- 
tion is a regular lattice, and derive the small k behavior 
of the SF of the daughter distribution. In Sect.|V]we then 
derive our result for the more general case that there is 
also arbitrary correlation between the displacements of 
daughters of different mothers, and then also consider 
the specific case of an initial regular lattice. As an exam- 
ple we derive the small k behavior of the SF of a lattice 
of correlated dipoles. In the final section we summarize 
our results and discuss briefly both possible developments 
and further applications of the results reported here. 



II. STRUCTURE FACTOR (SF) OF A 
STOCHASTIC CLOUD PROCESS 



We start with a spatial distribution of M "mother" 
particles in a cubic volume V, for which we write the 
microscopic particle density as: 






FIG. 1: The figure represents pictorially the generation of a 
a "cloud process" , as described in the text and characterized 
by Eq. (U). Starting from a "mother" point process (white 
particles) , a new particle distribution is generated by the "ex- 
plosion" of each mother particle into m (here, m = 3) daugh- 
ters each of which is displaced from the original point by a 
stochastic displacement u. The problem we solve is the de- 
termination of the two-point correlation properties of the new 
particle distribution given those of the "mother" distribution 
and the statistical properties of the displacement fields. 
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Let us denote by (...) the average over the ensemble of 
realizations of this point process, which we assume to 
be uniform at large scales, i.e., to have a well defined 
(positive) mean particle density ng in the limit 1^ ^ oo. 
The SF (or power spectrum) is then defined (see, e.g., 
i) as 
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where the limit V -^ oo is taken at fixed uq, and 
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is the Fourier transform (FT) of the point distribu- 
tion in the volume V. Note that, defined in this way, 
5„(k ^ oo) = 1, which is the usual normalization of the 
shot noise present in any particle distribution at short 
distances^. We suppose further that the point distribu- 
tion is statistically spatially homogeneous. In this case 
the SF is related by a FT to the usual two point corre- 
lation function. More specifically, if we denote by /i„(x) 



the off-diagonal part of the reduced two point correlation 
function, we have 

(7i(xo)n(xo + x)) -nl = no(5(x) + nlhn{x.) 

'^'^ S,,{k)e'^--. (3) 



no 



In other words 
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S'„(k) = 1 + nohnik) 

with hn{k) ^ FT[hn{x)]. 

Each particle in this distribution is a "mother" point 
which splits into to > 1 particles ( "daughters" ) . We take 
these latter to have identical unitary mass ■^. The daugh- 
ter particle distribution, which we call a cloud process, 
thus clearly has an average particle density po = rrniQ. 

Let us denote by Ur the displacement from the mother 
position of the r*'' particle (1 < r < to) belonging to the 
cloud generated by the «*'' mother (1 < i < M) of the 
distribution n(x). The resulting microscopic density of 
particles is then (see Fig. [1]) 



m M 
r— 1 z— 1 



(4) 



^ An alternative quite widely used normalisation (e.j 
mology literature) differs by a factor of I/tiq. 



^ We could equally take each daughter to have a mass 1/m, so 
in the cos- that mass is explicitly conserved. We choose unitary mass as 

this overall normalization factor cancels out in the SF. 



Let us now take the FT of p(x) in the volume V 
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We can now write the SF for /3(x) as 



5p(k) = hm 



|p(k;F)P 
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where (...) is, as above, the average over the reahzations 
of n(x) and (...) is the average over the reahzations of 
the displacement field. We will always assume here that 
the displacement field is statistically independent of the 
realization of the "mother" distribution, so that these 
two averages commute. In general, we can write 
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|p(k;F)P = Ar + ^^e-''-("^''-"^'') 
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where X^r^^s denotes the sum over r and s excluding the 
diagonal terms r — s, and similarly for X^Mi' 



III. UNCORRELATED CLOUDS: GENERAL 
RESULT 

In this section we assume the following statistical prop- 
erties for the displacement field: 

• The displacements applied to daughters with dif- 
ferent mothers are statistically independent; 

• The displacements applied to different daughters 
with the same mother may be arbitrarily corre- 
lated. 

More precisely, let us denote by p{u) the probabil- 
ity density function (PDF) for a single displacement, 
by ps(u, v) the joint PDF of the displacements u and 
V applied to two daughters of the same mother, and 
by pd(u, v;x) the joint PDF of the displacements ap- 
plied to two particles belonging to different mothers 
separated by x. Our assumptions imply that we have 
Pd(u, v;x) — p(u)p(v) for any x ^ 0, while we allow 
Ps(u,v) ^p{u)p{-v). 

Note that in writing the PDFs in this way, without 
labels for the clouds and for the particles in a single 
cloud to which the displacements apply, we assume im- 
plicitly the following symmetries: (i) p(u) is the same 
for all displacements, and (ii) Ps(u, v) does not depend 
on the cloud and is the same for all m{m — l)/2 cou- 
ples of the m particles belonging to the same cloud. In 
other words, if call P(ui,...,Um) is the joint PDF of 
the displacements applied to the m daughter particles 



in a cloud, we assume that P(ui, ...,Uj„) is the same 
for all clouds and is invariant under any permutation 
of the m displacement variables. This implies in par- 
ticular that Ps(u,v) — ps(v, u) and, as we show below, 
important constraints on the displacement-displacement 
correlations inside a single cloud. We clearly also have 
the consistency condition 



d'^vpsiu,v) =p(u). 



With these assumptions it is simple to show that it 
follows from Eq. ^ that 



|p(k)|2 = N + N{m-l)ps{k,~k) 
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where p(k) = FT[p{u.)] (i.e. the characteristic function 
of the single stochastic displacement) and ps(k, — k) is 
the following diagonal two-point FT of ps (u, v) : 



p,(k,-k)= fd''ud''vps{u,^)t 



-?'k-(u— v) 



Note that 

PsiK -k) - 0.(k) = y d''w;0,(w)e-*-- , 

where 0s (w) is the PDF of n-th relative displacement 
w = (u — v) between two particles in the same cloud: 



(jjsi'w) ~ / d ud ups(u, v)(5(w - u + v) . 

Since ps (u, v) = ps (v, u) , 0s (w) is an even function of w, 
and, as a consequence, 0s (k) is a real function. Further, 
from the fact that 0s (w) is a PDF, it follows that 0s (k = 
0) = 1, and |0s(k)| < 1 at all k. 

To perform the average (...) (over the realizations of 
the mother distribution), we use that 
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and that from Eq. ^, for asymptotically large M (i.e. 
for y — > c» with no fixed), one has 



/1,M 



Y. e-'"'-^^'-^^^ ) = M[5„(k) + (2^)'^no<5(k)] 



Using these results, together with N = niM and po = 
mnQ, we finally obtain 

Sp{k) = 1 + (m - l)0s(k) - m|p(k)|2 + m|p(k)|25„(k) . 

(8) 



This is the principle resuh of this section. It is useful to 
write it also in the form 



where 



5p(k) = 5(0)(k)+m|p(k)|25„(k). (9) 



5W(k) = 1 + (to - l)0,(k) - TO|p(k)|2 (10) 



depends only on the statistical properties of the displace- 
ment fields (i.e. independent of those of the mother dis- 
tribution). 

We note that for the case m — I Eq. ([5]) reduces to 



5i(k) = l-|p(k)|2 + |p(k)|25o(k), 



(11) 



This is precisely the equation derived in [i|, [^ for the 
transformation of the PS of a point particle distribution 
from S'o(k) to 5*1 (k) when each particle is randomly dis- 
placed, independently of the others, with a PDF p(u) for 
the single random displacements. 

Another simple case is that in which the mother par- 
ticle distribution 7i(x) is itself completely uncorrelated 
with hni'x.) = 0, i.e., generated by a homogeneous Pois- 
son process. Then S'„(k) = 1 and consequently Eq. ([S]) 
gives 

5,(k) = l + (TO-l)^,(k), 

Since ^s(k ^ 0) = 1 it follows that Sp(k -^ 0) = to, 
i.e., at large scales the cloud process is identical to the 
original distribution (up to a change in the mean parti- 
cle density by the factor in). This simply translates the 
fact that if a point process has no correlation, we cannot 
create correlation at large scales by dividing the particles 
into clouds by a stochastic process which incorporates no 
correlation between the different clouds. 



A. Properties of S'p(k) 

By definition any SF, and therefore the one we have 
derived, must satisfy the conditions 



Sp(k) > 
Sp(k -^ oo) = 1 , 



(12) 
(13) 



These must hold for any input 5„ (k) (itself obeying these 
properties) and any finite value of to,. 

The second property Eq. (|13p is simple to verify. It 
follows from the fact that both (/'^(k) and |p(k)| vanish 
in the large k limit. This is the case because they are 
FTs of functions which are integrable at the origin. 

That the first property Eq. p^ must be satisfied by 
our result is trivial: we obtained 5'p(k) by simply aver- 
aging Eq. ([6]) which is non-negative definite by construc- 
tion. However, as we now discuss, it is not simple as one 
might anticipate to verify it directly from Eq. ([5]). The 
property Eq. (J12p in fact encodes in a concise and very 



non-trivial manner constraints on the joint PDF Ps{u, v) 
which follow from the assumption that it is unique for 
any two particles in a cloud. 

Firstly we note that Eq. IT^ holds in fact if and only 
if 



5f(k)>0 



(14) 



for all k. This is the case because Eq. P^ must be 
true for an arbitrary S'„(k), and it is always possible to 
choose a mother point process for which it vanishes at any 
given k (taking, e.g., an appropriate regular lattice). For 
TO = 1 this condition is trivially satified, as |p(k)| < 1 
by definition (as FT of a PDF). For ?7i > 2 it may be 
rewritten as the condition 



M'^) > - 



1 



771—1 



-j\pi^)\' 



> -1 



(15) 



This does not trivially follow from the fact that p(k) 
and 0s (k) are FT of PDFs. As noted above, the latter 
gives only the weaker condition (/is(k) > —1. Clearly 
Eq. (|15[) encodes a non-trivial constraint on 0^ (k) , arising 
from the fact that it is related to a PDF for the joint 
displacement PDF. The latter is not simply an arbitrary 
normalizable function. The condition Eq. (jlSp tells us 
that we have in fact constrained it mathematically by 
the assumption about it we have made in our derivation: 
we have assumed that the cloud is generated in such a 
way that the joint two-displacement PDF is identical for 
all couples of particles. 

To illustrate this more explicitly we derive now the 
form taken by the constraint, in the form of Eq. ^ 

when Sp (k) is expanded in Taylor series around k 
Such an expansion can be made assuming that both p(u) 
and 05 (w) are rapidly decreasing at large u and w so that 
their FT are analytic at k = 0. We then have 



0. 



p(k)=^(-*) 



;(k- 



(16) 



1=0 



and 



,(k) ^p,(k,-k) = j^HY^^Ti^ ■ (17) 



;=o 
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The condition (fH|) can therefore be rewritten as 

4"'(k)- E^{("^-l)[k-("-v)]' (18) 

1=1 

-m ^ n liTTo^ X oT^I > , 

('=0 ^ ^ J 

where we have written explicitly w = (u — v) and used 
the symmetry assumption (k • v)-' = (k • u)-? for any j. 
This inequality, valid for all k, fixes all the constraints on 
the two-displacement correlation function in any cloud. 



Using again the fact that <?!>s(w) = 0s(— w), and mak- 
ing the further assumption that p(u) = p{—u), aU the 
odd power terms in this expression vanish so that we 
obtain: 



4"^(k) = ^^|(m-l)[k.(u-v) 



121 



(19) 
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E (2', )(k •")'<'-'') x{k-vr|>0. 



;'=o 

The leading term dominates at sufficiently small fc and 
therefore has to be non-negative. This implies the fol- 
lowing constraint on the correlations of any two displace- 
ments in the set of m correlated random displacements 
in each cloud: the matrix 



(m - l)-u(^)i;('^) + «('')«('') , 



(20) 



where \x,v ~ \^ ..., d, has to be non-negative definite. In 
our analysis in the next section, of the case that the 
mother distribution is a lattice, we will see how this con- 
straint, and ones which can derived at subsequent order 
in this expansion, simplify and are explicitly verified in 
certain cases. 



B. Behaviour at small k 

Let us now consider specifically the properties of the 
SF of the mother and daughter distributions as k — > 0. 

We note first that, because of the normalization condi- 
tions on the PDFs of the displacements, both 0s (k) and 
|p(k)| converge continuously to unity as k — > 0. It follows 

from Eq. ^ that S'^"^(k ^ 0) - fc", where a > 0. Sup- 
posing now that the initial (mother) point distribution 
has S'„(fc ^ 0) '^ fc''', we can infer that Sp{k ^ 0) ^ k'^ 
where (i) 7' = 7 for 7 < 0, and (ii) 7' = min{'y, a} for 
7>0. 

Thus the exponent of the SF around fc = can never 
be larger in the cloud process than in the original mother 
point process. Further it may differ from it (in which 
case it is smaller) only if S'„(fc = 0) = 0. Note that these 
conclusions hold independently of any assumption about 
the cloud process, other than that there is no correlation 
between the displacement sets creating different clouds 
and that the displacements are symmetrically distributed 
as shown above. 

This result may be explained more physically as fol- 
lows. The exponent of the small k behavior of the SF 
can be considered as a measure of the degree of order 
in the stochastic point process at asymptotically large 
scales [l5|, [l3l . The greater the exponent the more or- 
dered is the distribution. Indeed any lattice, which is the 
class of the most ordered particle distributions, the SF 
vanishes identically around k = 0, i.e., we can consider 
it to correspond to the behaviour k°°. Clearly a cloud 
process, without any correlation between the arrange- 
ment of matter in the different clouds, cannot increase 



the degree of order. That it may, on the other hand, 
decrease the degree of order when the mother distribu- 
tion has the property that 5„(fc = 0) = 0, reflects the 
difference between this class of distributions and those 
with Sn{k = 0) > 0. Indeed this difference is that under- 
lined by the classification of the former point processes 
as super-homogeneous (or hyper-uniform): the rapid de- 
cay of the density fluctuations at long wavelengths which 
characterizes them are the result of a delicate balance 
between small scale and large scale correlations in direct 
space. Indeed the condition ^^(fe = 0) = is explic- 
itly an integral constraint on the two point correlation 
function over all space. The processes which we are con- 
sidering, in which each particle "explodes" independently 
of all others, can break these global constraints by modi- 
fying only the small scale correlation properties. Instead 
for distributions with Sn{k = 0) 7^ 0^, which do not 
present such a correlation balance, such an uncorrelated 
re-distribution of matter at small scales cannot modify 
the nature of the system at large scales and also. 



Let us now consider the problem of the construction of 
a point process with a target behavior of S'(k -^ 0). From 
the results we have just derived it follows that a cloud 
process of the type we have just analyzed, without corre- 
lations between the displacements of the members of dif- 
ferent clouds, may be used to generate a distribution with 
the target exponent a > provided we start with 7 > a. 
In practice we can start with 7 = 00 by taking a regular 
lattice, for which 5n(k) = identically in a flnite region 
around k — (specifically, in the first Brillouin zone). 
The generated distribution will then have the exponent 
a, which depends through Eq. (jTU]) only on the statisti- 
cal properties of the displacement fields. Since a > the 
generated process is necessarily super homogeneous. The 
question we now address is what values of a are attain- 
able, and for what conditions on the number of daughters 
m and the displacement fields they are realized. We note 
that in [l| it has been shown that, for m = 1, one can 
realize by appropriate choice of p(u) any exponent in the 
range < a < 2. The upper bound a = 2 results for any 
p(u) with a finite variance, while the lower exponents are 
realized for PDFs with appropriately divergent moments 
of displacements at order less than two. For the case of 
correlated displacements, with Gaussian statistics, it has 
been shown also in [l[ (see also [1^) that a maximal value 
of a = 4 may be attained. We will see now that, with an 
appropriate value of m and conditions on the displace- 
ment fields, arbitrarily large positive target value of a are 
attainable. We will show that to obtain a certain value 
of a requires that one fixes a sufficiently large number of 
mass moments of all clouds of particles with respect to 
their respective initial position on a regular lattice. 



* These can be classified into Poisson-Hkc for Sn{k = 0) < -l-oo, 
and long-range correlated for Sn{k — > 0) ^ +00. 



C. Explicit expansion around fc = 

For this study of the smaU k behavior of >5'p(k) we 
consider the expansion of both (ps (k) and j'(k) in power 
series of k, which are given respectively by Eqs. (fT7|) and 
(fTH]) . Using these expressions in Eq. ^ we obtain 



1=0 

E°° , .,,(k • u)' 



[k.(u-v)]' 



1=0 



l\ 



l\ 



[5„(k)-l]. (21) 



Ah these expressions are vahd only under the assump- 
tion that all the moments in the sums are finite, while 
the derivation of Eqs. ([5]) and pHl) required only the in- 
tegrability of the probability distributions. If the prob- 
ability distributions p(u) and (/)s('w) have only a finite 
number of finite moments, the corresponding sums must 
be terminated at the appropriate order. There is then an 
additional term, of which the leading singular part can 
easily be determined. We will not discuss here the case in 
which there are such singularities ^ . The required gener- 
alization of the analysis described here is straightforward 
following the procedure defined in [l|. 



IV. UNCORRELATED CLOUD LATTICE 

Given the above motivation we now analyze in detail 
the cloud processes in the previous section for the partic- 
ular case that the mother particle distribution is a regular 
lattice, i.e. 



i(x)=^J(x-R) 



R 



where R is the generic lattice site. In this case the SF of 

n(x) is 



5„(k) = (2^)'^no^,5(k-H) 



(22) 



where the sum is over all the vectors H of the reciprocal 
lattice but H = 0. Note that this vanishes identically 
in the first Brillouin zone, and therefore in this region of 
the k— space the following relation holds exactly: 



5,(k)^l+(m-l)^(-0'(^ 



1=0 



E(- 

1=0 



:(k-u)' 



(23) 



^ We will see below that to obtain exponents greater than four 
with uncorrelated clouds the PDFs for the displacements must, 
in fact, have compact support. 



To simplify the presentation of our determination of the 
conditions required to have an arbitrary (analytic) small 
k behavior of S'p(k), we start with the one-dimensional 
case, for which Eq. ([^5]) becomes 



Sp{k) 



-1) > (-i/c)'^ — TT-^-m 



II 



E(-*^)' 



II 

(24) 

where u and v are the displacements applied to two differ- 
ent particles belonging to the same cloud. First of all we 
see immediately, as above underlined, that the final par- 
ticle distribution p{x) is superhomogeneous, as the zero 
order contribution in k to Sp{k) vanishes identically for 
any choice oi Ps{u,v) [i.e. of 4>s{w)]- In the notation of 
the previous section, we have explicitly that a > 1. We 
can write 



i:<-'^)'t 



;=o 



E 

/=0 



i-^k) 



n 



I I 

E(-i: 






where we have used u and v instead of only u as the mo- 
ments of the single displacement are the same for every 
particle. Moreover, by expanding the terms {u — vf in 
Eq. ((241), we have 



E(-^^)' 



{u — vY 



1=0 

Therefore Eq. 



II 



OO / ., N I I 



/=0 

becomes 



3=0 



ui X w'"-? 



(-^fc)' 



1 = 1 ■ j=0 



(to — l)wJ X v^~i — mui x t)'~J 



(25) 



Making the additional assumption of statistical symme- 
try in the displacements, p{u) = p{—u), all the terms 
with odd / in Eq. (|25p vanish. 



A. Order by order analysis and conservation of 
mass moments (d = 1) 

Let us now analyze in detail Eq. (|25p . denoting by 
0„(fc) its term proportional to fc". Given our hypotheses 
the lowest order non-zero term is n = 2: 



02 (fc) = V? + (m — l)u X V k 



It 



is simple 



to 



verify explicitly that 



u^ + (to — 1)u X V > always, as required (from 

the fact that Sp{k) is a SF). First of all it is the 
one-dimensional version of the condition ([20|) . This can 
be seen more directly as follows: if we denote by Ui with 
i = 1, ..., TO the displacements applied respectively to the 



m daughter particles originating from the same mother 
(i.e. belonging to the same cloud), it is clear that 



E^ 



> 0. 



This quantity, however, given our symmetry hypotheses 
about the displacments distribution is nothing other than 




+ (m — l)u X V 



Consequently the condition to have an identically van- 
ishing 02 (fc) term, and therefore a small k SF of order 

greater than two (i.e. a > 2), is (X^I^Li ""0 = ^^ or in 
other words, 



E^ 



= 



(26) 



with probability one^. This means that the center of 
mass of each cloud does not move away from the mother 
particle when the displacements are applied. Clearly for 
771 = 1 this condition can only be trivially satisfied by ap- 
plying no displacement, in which case the daughter dis- 
tribution is the original lattice distribution. For m — 2 it 
can be satisfied non-trivially: choosing the displacement 
of a first point with the PDF p{u), the other particle 
is then displaced deterministically by —u. For m > 2 
the condition can be satisfied while admitting a higher 
degree of stochasticity: it fixes deterministically the dis- 
placement of only one particle among m once the other 
(m — 1) are chosen stochastically. 

The analysis of the term 04{k) is more complex. We 
will now show that the condition for it to vanish is one on 
the second moment of the mass dispersion of each cloud. 
Directly from Eq. pS)) we have 



Oi{k) 



4 r 



k 
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_ (27) ■ 

where we have denoted a^ — u^ and again used the as- 
sumed symmetries of the displacement field. If the term 
of 02{k) vanishes at all k, i.e., if Eq. (P5|) is satisfied for 
each cloud, then the O^^k) term must be non- negative 
[since the SF Sp{k) must be non-negative at all fc]. Thus 
the coefficient of fc^ in Eq. ((27|) must be non-negative. In 
order to show that this is indeed the case we note first 
that Eq. ((26|) implies that 



Ui 2. ""i 



- „4 



-|- (777 — 1)^-^1) = 



^ Note that if we had allowed an asymmetric p{u) with a non-zero 
average value u the same condition 02{k) = would have been 



By using this relation in the coefficient of Eq. 
then have 



23 we 



—u^ — 4(777 — l)u^v + 3(777 — 1)77^7;^ — "ima 

= 3w^ + 3(777 — 1)77^U^ — 3777CT"' . 

Further it is simple to show that 

3777(77* + (777 — l)77-^7;-^ — 777CT ] 



(28) 



T.^] - E 



u' 



which is manifestly non-negative. 

This result also gives the condition necessary in order 
to have both the 02{k) and 04{k) identically vanishing: 
once Eq. ([26]) for the conservation of the center of mass 
of each cloud of the system is assumed, in order to make 
the variance in Eq. (j28p vanish one requires also that 



(29) 



E^ 



for each cloud, with probability one. 

In summary: in order to obtain with this algorithm a 
particle distribution with a SF of order larger than fc^ 
at small k, one has to satisfy exactly the following two 
conditions: 

• every cloud must have the same displacement of its 
center of mass from the initial point of the cloud or 
mother position [and, in particular, iip{u) = p{—u) 
the center of mass of the cloud must coincide with 
the position of the mother particle] ; 

• every cloud must have the same inertial moment 
(or second moment of its mass dispersion) with re- 
spect to the initial point. The value of this inertial 
moment is fixed by the second moment oi p{u) and 
m as ma^ . 

This analysis can be continued in order to determine the 
conditions needed to obtain a Sp{k) of order higher than 
fc^" at small fc, for any integer 77. The result is simply 
that in order to obtain this goal one has to fix the first 77 
moments of the mass dispersion 



E-H 



with j = 1 , 



(30) 



written as JZiS 



: mu for all the clouds. 



where the constants Cj are determined by the j-th mo- 
ment of p{u) and 777 as Cj = mu^ . Clearly this gives 
n conditions and consequently one has to have at least 
m = n particles in each cloud in order to make the re- 
quirements given by Eq. pop realizable. For m = n there 
may be a single non-trivial solution to the constraints, 
i.e., a unique choice of displacements. In this case the 
generated distribution will be a lattice with basis, with a 
SF which again vanishes in a finite region around fc = 0. 



For m > n the set of constraints may be satisfied (for 
some range of values of the constants Cj) while leaving 
free (to — n) degrees of freedom. These may be then fixed 
stochastically, leading generally to a stochastic particle 
distribution ^ with a leading non-zero term at C2(n+i)(fc)- 
For n = 2, for instance, we need at least two particles. 
In order to fix the center of mass of the pair at the lat- 
tice site, and its second mass moment to a given value C2, 
clearly fixes deterministically the points to lie at iycT/S. 
Taking three particles one can instead satisfy the con- 
straints fixing one degree of freedom stochastically: plac- 
ing one point at u with probability p{u), the position u' 
of a second point is determined by solving the quadratic 
equation (m')^ + uu' + (u^ — C2/2) = 0, and finally a third 
point is placed at —{u + u'). Note that the existence of 
a solution to the quadratic equation places a strict up- 
per bound on M, M < ■\/2c2/3. Thus the probability p(u) 
necessarily has finite support (and cannot in particular 
be Gaussian) which is proportional to y^. It is clear 
that this is a general requirement for any algorithm of 
this kind producing a SF with a leading small k behavior 
fc" with n > 4: in order to make the coefficient of the k^ 
term vanish the second moment of the mass dispersion 
of the cloud must be limited with probability one. Since 
the displacement of any particle contributes in propor- 
tion to its square, the probability of displacement larger 
than JII2 must be zero. 



B. Generalization to d > 1 

In dimensions higher than one the problem is essen- 
tially the same. The analysis is, however, considerably 
more complicated because of the vectorial nature of the 
displacements. 

Let us first consider the conditions required to make 
the terms 0{k^) and 0{k'^) vanish. Fixing the center of 
mass in d dimensions gives d scalar equations 



i=l 



.ip) 



where u^ is the fi (with fx — 1, ...,d) component of the 
displacement of the i*'* particle of the cloud. To satisfy 
this condition non-trivially evidently requires that there 
are at least two particles in each cloud. Fixing the second 
moments of the mass dispersion of the cloud gives d(d + 
l)/2 scalar equations, i.e., d equations of the form 



E 

i=l 



,(m) 



^fj.fi 



> for y^ = 1, 



'' Here we mean by "stochastic" that there is a non-empty compact 
domain of k— space in which S'p(k) is continuous and strictly 
positive. 



with a^^ ~ TO,[u^^ ]2, and d{d — l)/2 equations of the 
form 



y^ uf'u^ ^ a^,i, for 1 < yu < J/ < d , 



with a^^ = ?Tiu(^)-u(''). Therefore to obtain a SF 5p(k) 
of order larger than the fourth at small k imposes 
[d + d{d + l)/2] = d{d + 3)/2 scalar constraints on the 
displacements. This counting of constraints may be con- 
tinued to higher orders, determining the number of con- 
ditions N{n, d) which must be satisfied to obtain, in d 
dimensions, an S'p(k) vanishing faster than fc^" at small 
k. Noting that all the moments of given order, say /, of 
the mass dispersion constitute a fully symmetric Z-rank 
tensor in d dimensions, which has ( "*" ~ ) independent 
components, we find 



N{n,d)^Y. 
1=1 



d + l-1 
I 



Generalizing the reasoning for the case d — 1, one might 
then be tempted to conclude that, since each particle 
brings d degrees of freedom, the minimal number m of 
particles per cloud required to given a 5'p(k) of order 
larger than /c^" at small k is M{n, d)/d. This conclusion 
is, however, not correct in that it tells us the number of 
particles required to satisfy such conditions for arbitrary 
physical values of the mass moments. To make the coeffi- 
cients in the expansion of S'p(k) vanish up to some order, 
while remaining non-vanishing at subsequent orders, we 
require only that a number of particles sufficient to al- 
low us fix a set of physical values of the mass moments 
up to a certain order, while allowing higher moments 
to vary. Put another way, by imposing additional sym- 
metries or constraints on the PDF of the displacements, 
we can reduce the number of non-trivial constraints (i.e. 
equations), reducing the others to simple identities. The 
following example illustrates this point trivially: in d > 1 
we can always make a cloud lattice by putting together 
at infinite number of one dimensional cloud lattices, i.e., 
by constraining the displacements of the particles to lie 
along a chosen axis of the lattice. The calculation in 
d = I then remains valid, as all the additional constraints 
on the mass moments of the clouds with components in 
orthogonal directions are trivially satisfied. The same is 
true in fact if the displacements are in an arbitrary (but 
fixed) direction. The number m of particles per cloud 
required to obtain a SF with given leading order then 
remains the same as in cf = 1. One can also evidently 
consider less radical "dimensional reductions" , taking in 
d dimension the displacements of particles in each cloud 
only in a hyperplane of dimension smaller than d. 

Even without such a reduction to a lower dimensional 
problem, it is easy to give examples in d > 1 which sat- 
isfy the constraints required to make all terms of the SF 
up to a certain order 2n vanish, with much less than 
M{n, d)/d particles, and which do not have the feature 
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FIG. 2: The figure represents pictorially the generation of 
a a cioud process starting from a lattice, in d = 2 with a 
"dimensional reduction" to d = 1 in which the two particles 
in each cloud are displaced along the same direction. 



FIG. 3: The figure represents pictorially the generation of a 
a lattice "cloud process", in d = 2, in which each cloud is 
symmetric with respect to its lattice site. 



CORRELATED CLOUDS 



of the above examples of breaking the statistical spatial 
isotropy of the clouds. Consider, for example, the case 
that each cloud contains an even number of particles, ar- 
ranged symmetrically with respect to the center of mass. 
From the derivation we have given above, it follows im- 
mediately that SpCk) (X k^ at small k with 7 > 4 (as the 
inversion symmetry gives automatically the conservation 
of the centre of mass) . As shown above the coefficients of 
the term proportional to k^ will vanish if the second order 
moments of the mass distribution are cloud- independent. 
This can be attained, for example, in d dimensions, by 
taking d pairs of particles arranged symmetrically about, 
and equidistant, from the origin, and all mutually or- 
thogonal. This gives a second moment of the mass dis- 
persion of the cloud which is proportional to the identity 
matrix, and therefore invariant under a random rotation 
R £ SO{d) of the whole configuration . Further it is pos- 
sible to show (see Appendix A for detail) that, because of 
the imposed inversion symmetry, the terms proportional 
to k^ vanish identically. Thus, placing such a cloud with 
an orientation chosen randomly at each site, one obtains 
a leading non-zero term of order k^ . This term is non-zero 
because the (tensor) fourth mass moment is not invari- 
ant under rotation of the configuration. Further, if the 
stochastic process determining the orientation is statis- 
tically isotropic, the SF at small k reflects this isotropy 
and is a function of k only (rather than the vector k) . It 
is simple to check that the number of required particles 
for this algorithm 2d is less, for any d, than the number 
A/'(3, d)/d given above. 



We now consider the case in which also displacements 
applied to particles belonging to different clouds may be 
arbitrarily correlated. In the notation introduced above 
for the joint PDFs of the displacements, this means that 
we now assume, at least for some separation of mother 
particles x, that Pd(u, v, x) ^ p(u)p(v). Further we make 
the following natural assumptions: 

Pd(u,v;x) =pd(v,u;-x) 
lim pd(u, v; x) = p(u)p(v) . (31) 

cc— *oo 

We have, of course, also the consistency condition 



d'^fPd(u,v;x)-p(u) Vx (32) 



Note also that, strictly speaking, the displacement field 
is not continuous as a function of spatial separation be- 
cause the correlation u • v(x) between two displacements 
applied to two spatial point at vector distance x in gen- 
eral does not converge to u^ for x ^ 0. The natural 
choice for this limit is Pc((u, v;x ^ 0) = ps(u,v). How 
precisely this limit is taken will not, however, be of im- 
portance for our main application concerning a regular 
lattice distribution of the centers of the clouds, as in this 
case the distance between two different centers is always 
different from zero. 

In order to find the expression for S'p(k), once S'„(k) 
and Pd(u, v; x) are given, we have to go back to Eq. ^. 
From this it is simple to show that 



ip(k)r 



N 



l + (m-l)p^(k,-k) 

2 ^'^ 
™ ^^-,k.(x.-x,)p^(k,-k;x,,), 



^¥"3 



where x^ = x^ — x^ and 

PdiK -k; x) = / d'^u d'^ypdiu, v; x)e-^'^("-^) . 

Note that, analagously to ps(k, — k), the function 
"PdO^, — k; x) is the characteristic function of the stochas- 
tic vector w(x) = u(x) — u(0), i.e., of the difference 
between the displacements apphed to two particles be- 
longing to clouds whose mother particles are separated 
by X. In other words, if '/'d(w; x) is the PDF of w(x), we 
have 

Pd(k,-k;x) = 0d(k;x), 

where (/)c;(k;x) = FT[(pii{w;x)]. In order to perform the 
average (...), we recall that in general for any function 
/(x, y) of two spatial variables one has 



/1,M 

E 



/(x.,x,)\-n2yy^d'^a;A/(x,y)[l + /i„(x-y)] 



where /i„(x) is the two point correlation function as de- 
fined in Eq. ^. Using this relation, we obtain 



/1,M 

E' 



-'ik-(xi— Xj) 



0d(k;x,j) ) =nl d'^xd'^y x 



,— 'k-(x-y) 



i(k;x-y) [l + /i„(x-y)] 



(33) 



Taking the limit T^ -^ oo with N/V = po and m > 2 
fixed, we arrive at our result: 

5p(k) = l + (m-l)0,(k) (34) 

+P0 / d'^x e-'^'-^M^; x) [1 + /^„(x)] - {27rfpoSik) . 



Note that the case m = I can be included in Eq. ((34|) 
by considering, only for this value of m, Pd(u, v; x) to be 
spatially continuous, i.e., converging to p(u)(5(u — v) for 
X ^ 0. The result Eq. ([M)) for m — 1 then agrees exactly 
with the analagous equation found in [l| for the PS of a 
particle distribution after the application of a correlated 
displacement field. 

We will not analyse here the small k expansion of 
Eq. (j34p . but simply note that such an analysis can be 
done easily by following the steps for the study of the 
similar equation for the case m = 1 in p|. 

Let us now give the special expression of Eq. (|M|) for 
the specific case in which n(x) is a regular lattice point 
distribution. 



Correlated cloud lattice 
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where R runs over all the lattice vectors except R = 0. 
We therefore have 

Sp{k) = l + (m-l)<^,(k) (35) 

+ mJ2 e"*-''<^rf(k; R) - (27r)'^po^(k) . 
Rt^o 

This formula is a good starting point for the study of the 
small k behavior of Sp{k) for a cloud lattice for different 
choices of pd(u, v;R) and m. We now use the following 
chain of identities to rewrite Eq. ([34]) in a more useful 
form: 



i2n)''poSik) = (2^)'^po 



Y,6ik-n)-J2d{k~li) 



H 



^E' 



R 



JkR 



mSn{k) , 



(36) 



where we have used the definition ((22|) of S'„(k) for a 
regular lattice and the lattice identity 



Y^ e-*-^ = {2TTyno J2 '^(k - H) , 

R H 

By using Eq. ([36]) . we rewrite Eq. ([35]) as 



(37) 



Sp{k) = (m-l)[0,(k) 



1 



(38) 



-I- m y^ e 

R^^O 



ik-R 



(k;R)-l +jnSn{k) 



Note that 5„ (k) vanishes identically in the first Brillouin 
zone, so that it does not contribute to the small k expan- 
sion of Eq. (|55|) . Therefore, assuming all the moments of 
the displacements to be finite, in this region of fc— space 
we can write 



5,(k) = im-l)Y^{-^y^^:^ 



(39) 



^^1 

^E-^^-^'EH) 

Ti=iO 1=1 



[k-Wrf(R)]' 
l\ 



where we have denoted, respectively, w^. the relative dis- 
placement of two different particles belonging to the same 
cloud, and "Wc;(R) the relative displacement of two parti- 
cles belonging to two clouds whose centers are separated 
by the lattice vector R. From this formula we can de- 
duce the conditions on the two-displacement correlations 
to have a given small k SF for the resulting particle dis- 
tribution. 



When n(x) is a regular lattice 

/i„(x) = — V 5(x - R) - 1 

71 r^ ^ ^ 



no 



Rt^o 



B. Example: lattice of random correlated dipoles 

Let us consider the following example as an application 
of Eq. ([55]) . Each particle on a perfect lattice is split 
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into two particles, and the following displacements are 
applied: 



Ui(R) = +r/(R), U2(R) = -?7(R) 



(40) 



where r/(R) is a random vector at each R specified by a 
lattice translationally invariant correlated stochastic pro- 
cess. The average over all the realizations of the displace- 
ments of a function ^({^^(R)}) of the displacements may 
be written as the functional integral 

X = n / A(R)^ (Mm) ^({^(R)}) • (41) 

R "^ 

We assume that the joint probability density function of 
all the displacements V ({^^(R)}) is invariant under any 
lattice translation. Moreover it is simple to show, given 
our symmetry assumption for the displacements, that we 
can take V ({?y(R)}) to be invariant under the change of 
sign of any individual f](R). This ensures that Pd(u, v; x) 
is well-defined as required in our derivation, i.e., the joint 
probability for displacements to two particles at different 
sites is the same for all couples, and depends only on 
their relative separation. With this assumption it follows 
that all odd powers i in the sums in Eq. ((39)) vanish. 
Calculating the contribution at second order in k we have 



[k-w,]2^4[k-r;]2 



(42) 



and 



[k-Wd(R)]2 = yd3,/(0)d3r,(R)prf[,^(0),r;(R);R] 
x[k.(r;(0)-r;(R))]2 



2[k-r?]^ 



(43) 



The latter resul ts uses the fact that the two point corre- 
lation function rj{0) ■ t?(R) ~ 0, because of the assumed 
inversion symmetry. Using these results in Eq. (|39p to- 
gether with the identity Eq. [371 one finds that the two 
contributions cancel in the first Brillouin zone to give 
zero. 

At next non-trivial order, fourth order in k, we find 



[k • Ws]-* = 16[k • 7/]4 



(44) 



and 



[k.Wd(R)]4 = yd''77(0)d'^77(R)pd[77(0),r;(R);R] 

x[k • (7 7(0)-r;(R))f 
= 2[k • 77]4 + 6[k • r;(0)]2[k • ?7(R)]2 .(45) 

Using again the identity Eq. ([37]) . we obtain the leading 
non-trivial contribution to the PS which may be written 

Sp{k) ~ ^k^kpk^ks Y, e-''"^r;„ (0)7;^ (0)77^ (R)?7A-(R) , 

R 

(46) 



where we adopted the sum convention on the repeated 
index a, (3, 7, 6. The sum over R is just the FT (on the 
lattice) of a two point correlation function, the behav- 
ior of which as k ^ depends on the nature of these 
correlations. The leading small k behaviour of the SF 
will depend on that of the sum. If the correlations of the 
dipoles are short-range, the sum converges to a positive 
constant at small k, giving a leading behaviour propor- 
tional to fc^. If, on the other hand, they are long range 
correlated, this sum will diverge as a power-law at small 
k, with an exponent less than that of the dimension of 
the space. This will lead to behaviour of the overall SF 
proportional to k'^ with 4 — d < 7 < 4. Finally, if the 
correlations of the dipoles have themselves superhomoge- 
neous properties^, one can obtain such a behaviour with 
7 > 4. 



VI. SUMMARY AND DISCUSSION 

In this paper we have introduced and analyzed a wide 
class of non-trivial stochastic point processes for which it 
is possible to write exactly the two-point correlation func- 
tion and/or SF. They are obtained from a given "mother" 
particle distribution, the SF of which is assumed known, 
by substituting each particle with a cloud of a fixed num- 
ber 771 of other particles. The position of the new particles 
composing the clouds is determined by that of the re- 
lated mother particle plus a stochastic displacement vec- 
tor. An important assumption in all our calculations of 
the SF is that the stochastic process determining the dis- 
placements of the particles is independent of the mother 
distribution, i.e., the displacements of the cloud particles 
do not depend on the properties, statistical or otherwise, 
of the mother distribution to which they are applied. In 
practice this means that our SF is defined with respect 
to an ensemble average over two independent ensembles: 
one describing the realizations of the mother distribution, 
the other those of the displacement process. 

We have distinguished two cases in which: (i) only 
the displacements of different particles belonging to the 
same cloud may be correlated, and (ii) the displacements 
of particle belonging to different clouds may also be sta- 
tistically dependent. In both cases we obtain a direct 
generalization of the relations found in [IJ for m = 1. 
In the first case, once the average over realizations is 
taken, the SF of the final particle distribution is related 
to the SF of the mother distribution by a local relation 
in the wave vector k, while the second case leads to a 
more complex relation. A detailed analysis of case (i) 
led us to find and to discuss, in the case of an initial 
lattice mother distribution, the relations linking the ex- 



* Note that, as the functions ric,(Q)r]fj {0)ri-,f (Rjrjg (R) for a = 13 and 
7 = 5 are non-negative at all R,, their lattice Fourier transforms 
can vanish at k = only if they are identically zero for all R.. 
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ponent of the final SF to the number of conserved mass 
moments in each cloud: we have seen that such an ex- 
ponent, if larger or equal to 4, tells us directly which are 
the locally conserved mass moments in the distribution. 
When we move to case (ii), the presence of cloud-cloud 
displacement correlations "interact" with the local mass 
moments conservation in determining the small k behav- 
ior of the SF of the particle distribution. 

One application we have in mind of the results derived 
here for this class of "stochastically ordered" point pro- 
cesses is in the systematic study of the dynamics of parti- 
cle systems driven by long range pair interactions. More 
specifically, in the case of gravity, it is expected (see, e.g., 
[8|) that the large scale fluctuations in an infinite parti- 
cle system dominate the dynamics of the gravitational 
clustering for an initial SF with a small k behaviour pro- 
portional to fc^ and 7 < 4. This "hierarchical" behaviour 
has been observed numerically for a range of such 7, up 
to a maximal value of 7 = 2 (see [20] for a recent discus- 
sion, and further references) . No study of the regime of 
initial conditions 7 > 4 has been performed up to now, 
as no algorithm has been given in the literature, to our 
knowledge, which can generate such an initial condition^. 
It is expected that the gravitational clustering will be 
qualitatively different in this case, with structures being 
built up from smaller to larger scales. Indeed the rea- 
son why we expect such a difference can be understood 
easily in the context of our constructions here of such 
point processes. For 7 > 4 the fluctuations in mass are 
so suppressed that gravity is effectively "screened" (at 
least in the initial conditions): in a multipole expansion 
of the mass far away from a given point in the "uncor- 
related cloud lattice" only the leading moment varying 
from cloud to cloud will contribute (as the contributions 
from the moments which are fixed will cancel out). If 
this leading contributing moment is the second moment, 
this gives an effectively short-range force (decaying as the 
inverse of the fourth power of the distance) . 

Our analysis here has been limited to the case of "ana- 
lytic" exponents for the SF derived from short tailed dis- 
placements probability density functions (PDF). In prin- 
ciple one can consider also the case in which such PDFs 
have a long power law tail. In this case a singular part of 
the small k expansion of the related displacement char- 
acteristic function arises. This leads (see pj, i9j]) to a final 
SF characterized by non-analytic (e.g. fractional) small 
k exponents. We have seen, however, that in the uncor- 
related cloud lattice, to attain powers larger than 7 = 4 
we must in fact take limited displacements, and thus we 
necessarily obtain "analytic" exponents (and, in fact, a 
analytic behaviour of the SF at small k). In the case of 
the correlated cloud lattice, nevertheless, we have given 



an example (random correlated dipoles) which shows how 
such non-analytic powers should be attainable by includ- 
ing appropriate correlations between the clouds. 

We return finally to an important feature of Eq. ([M]) 
which we have discussed at some length in Sect. lIIIl This 
is that the exponent in the small k scaling behavior of 
5p(k) cannot be larger than that in S'„(k). As we ex- 
plained, this can be understood physically as it means 
that the replacements of particles by clouds cannot make 
the initial particle distribution more ordered. We under- 
line, however, that this is true given the assumptions we 
have made, and specifically assuming that the stochas- 
tic process generating the clouds is independent of the 
mother distribution. In a forthcoming article with an- 
other collaborator [23j we will report results on a re- 
lated kind of construction of superhomogeneous point 
processes, starting from tilings of space with equal vol- 
ume tiles. In this case it turns out that one can, in certain 
circumstances, ascribe a cloud of particles to each parti- 
cle of a given point distribution and as a result increase 
the exponent 7. The reason why this becomes possible 
is that the displacements of the particles are not applied 
independently of the correlation properties of the under- 
lying point process, as we have assumed here. Indeed 
in order to increase the exponent requires that the mo- 
ments of the clouds are "tuned" appropriately to the tile 
in which the mother point lies. We note also that, in the 
case that the initial tiling is taken to be a lattice, the 
algorithm described coincides with that given here and 
similar results to those given here are recovered. Details 
will be reported in [23| . 

We thank B. Jancovici, J. Lebowitz and S. Torquato 
for useful discussions. 



APPENDIX A: SYMMETRIC CLOUD LATTICES 

In this appendix we study the SF for an uncorrelated 
cloud process on a lattice, for the case that each cloud is 
symmetric with respect to its own lattice site. 

The number density of such a particle distribution can 
be written as: 



P(x) 






R 



.(R.) 



(Al) 



R 



where R is the lattice site and m is the number of parti- 
cle per cloud. Since we consider symmetric clouds, m is 
even and can be written as m — 2p where p is a positive 
integer. We consider initially a finite lattice with M sites 
(occupying a corresponding finite volume V), and then 
send M ^ CXI at the end of calculations. Taking the FT 
of Eq. (jAip , we have 



^ Explicit algorithms generating point distribution with 7 = 
4 have, on the other hand, been discussed. See, notably, 



liLaJUi 



p(k,y) = 5:5:e-'^(^+"r). (A2) 

R j=l 

Now we impose the inversion symmetry for the clouds 
with respect to their center of mass (i.e. the lattice site). 
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This means that for each particle at (R + u) there is 
another particle placed at (R — u). Therefore for each 
cloud we can count with j from 1 to p = to/2 a set 
of particles which are not the symmetric image one of 
each other, and with p + j their respective symmetric 
images. Particles in a single cloud with I < j < p can be 
arbitrarily correlated. Let us call as above Ps (u, v) the 
joint PDF of a couple of displacements referred to the set 
of particles with j = 1, ...,p in the same cloud and p(u) 
the PDF of a single displacement. 

Imposing this symmetry of the clouds we can rewrite 
Eq. 1X2]) as 



p(k, y ) = 2 ^ e-^"--^ ^ cos(k . uf ) ) . (A3) 



R 



Then taking the squared modulus we obtain 

\piKVr (A4) 



l,p 



4 ^ e-''^-^^-^'' ^ cos(k • uf )) cos(k • u 



(R') 



R,R' 






l,j 



^P)? 



i,p 
+ EE ^o^C' ■ ur^) cos(k • u|^^) 

R l^j 

+ J2 e-"^-(^-^') ^ cos(k • uf ) cos(k • u|^')) 

R^^R' 

We can now take the average over the displacements. 
In order to do this we recall that we are assuming that 
the displacements related to particles belonging to differ- 
ent clouds are uncorrelated. Consequently the first and 
the third terms in Eq. (jA4p have to be averaged over 
only p(u), while the second one as to be averaged over 
Ps(u, v) containing all the two-displacements correlators 
in a single cloud. This gives 

(|p(k,l/)|2)=4Mp([cos(k-u)]2) 
+4:Mp{p - 1) (cos(k • u) cos(k • v)) 
+4/ J2 e-^"^-'^-^'' (cos(k • u))' 

R/R' 

= 2iV ([cos(k • u)]2) + N{m - 2) (cos(k • u) cos(k • v)) 



-iVTO(cos(k-u))^ + Li. 



(A5) 



where N — Aim is the total number of particles and 
"l.t." indicates a "lattice term" which is proportional to 
the lattice SF (and which therefore does not contribute 
around k = 0). In performing the last step of Eq. (|A5[) 
we have used the simple identity 

y^ g-.k-(R-R') ^ y- g-Jk.(R-R') _ j^j 



which gives rise to the third term of the last step. Now 
we use the definition of the SF 



lim 

M^oo 



{\P{KV)\^) 



N 



which gives 

Spik) = 2([cos(k-u)]2) + (TO 



m (cos(k • u)) + l.t 



- {2nfpoSik) , 

2) (cos(k • u) cos(k • v)) 
(A6) 



It is simple to verify that, as expected, Sp{0) — 0. We 
now expand Eq. (|A6p in power series of k and study it 
order by order: 



\a+b 



Sp{k) 



^ (-1)' 

t^ (2a)!(25)! 

\2a 



2 ( (k • U 



,2(a+6) 



(A7) 



\2b\ 



-f (to - 2) ((k • U)2°(k • V)2'') - TO ((k • U)2'') ((k • U) 

We see that only even powers of k are present. Let us 
call 02n{k) the term of order fc^" in the series above. We 
see immediately that 

02(fc) = 

as all the clouds conserve the center of mass at their 
lattice site for symmetry. Therefore the first non trivial 
term is O^i^k), which after some manipulation can be 
written as 

0,ik) = i[((k.u)4)-f(p-l)((k.uf(k.v)2) 



- p{{^-nry 



(A8) 



Oe{k) 



[{{k-nf) + {p-l){{k-unk-^)') 



The next term is 
1 
"l2 

- p{{k-uf){{k-ur)] 

Both terms are in general non-zero. In order to find the 
condition for C'4(fc) to vanish, we rewrite it as 



(A9) 



04(fc) = i{([(k.u)2-((k.u)2)]^) (AlO) 

-F(p-l)([(k.u)2-((k.uf)][(k.vf-((k.vf)]}) 

In any cloud there are p = to/2 stochastic displacements 
which constitute a closed set of p symmetrically corre- 
lated variables (i.e. the correlation between any pair of 
these displacement is constant and there is no correlation 
with displacements in other clouds). Since the correlation 
matrix of the random variables [(k • Uj)^ — ((k • u)^)] , 
with j — 1, ...,p in a single cloud, has to be, as all the 
correlation matrix of a set of random variables, positive 
definite (see above) we always have 04 (fc) > 0. This can 
be seen in a more intuitive way by noticing what follows: 




J2[i^-n,r-{{k-u,f)] 



2pOi{k) (All) 
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which is a variance and consequently manifestly non- 
negative. Equation (|A11[) also implies that C'4(fc) van- 
ishes if and only if 



i.e., if with probability 1 we have 

f:(k.u,f =p((k.u)^). 

This is just the "conservation law" of the second moment 
of the mass dispersion of the clouds in the direction of k. 
As the orientation of k is arbitrary, this means that, in 
order to have 04 (fc) — identically, the second moment 
of the mass of the clouds must be conserved as a tensor 



/„ 



E? 



u]: U: With /i, ly 



We now analyze Oe{k) which is given by Eq. 
First of all we note that 



E 



(k-u,r 



p 

E 
1=1 



[(k' 



Ui 



((k-uO^)] 



-12pOe{k) . 



1, 



(A12) 
But, as seen above, if C'4(fc) = the second sum in 
Eq. (jA12[) vanishes identically, and therefore we can con- 
clude that when C4(fc) = identically also Oe{k) — 
automatically, and the dominant term in Sp{k) becomes 
Os{k). 

One can continue this analysis further and show, after 
some more involved algebra, that the dominant term in 
the small k expansion of the SF is of order fc^" with n 
integer, and next order terms are proportional to fc4n+29 
with q again an integer. The exponent n depends on the 
order to which the moments of the mass dispersion of the 
clouds are conserved. 
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